Nonequilibrium Rashba field driven domain wall motion in ferromagnetic nanowires 
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We study the efTects of spin-orbit interaction (SOI) on the current-induced motion of a magnetic 
(Bloch) domain wall in ultrathin ferromagnetic nanowires. The conspiracy of spin relaxation and 
SOI is shown to generate a novel strong nonequilibrium Rashba field, which is dominant even 
for moderate SOI. This field causes intricate spin precession and a transition from translatory to 
oscillatory wall dynamics with increasing SOI. We show that current pulses of different lengths can 
efficiently be used to control the domain wall motion. 
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Introduction. — The efficient and reliable manipulation 
of magnetic microstructures forms the basis for most in- 
formation storage devices used nowadays. Commonly, 
a local magnetic field controls the magnetization in a 
given domain. The alignment direction encodes a clas- 
sical bit, and the data storage density is limited by the 
domain size. Recent advances in nanofabrication have 
implemented Berger's proposal [H of moving a domain 
wall (DW) by a current-induced spin torque 0. This 
may allow for ultrasmall magnetic devices with high data 
storage density in itinerant ferromagnets carrying a spin- 
polarized current, where the local magnetization is elec- 
trically controlled by a spin torque arising from the ex- 
change coupling to conduction electron spins. Across 
a wide DW, the spin polarization adiabatically follows 
the magnetization direction, angular momentum is trans- 
ferred due to total spin conservation, and DW motion 
along the current direction is induced. In addition to 
this adiabatic spin torque 0] , a nonadiabatic spin torque 
occurs This "/3-term" is due to spin relaxation, 

which causes the itinerant spin polarization to "lag" be- 
hind the magnetization. Its contribution is essential for 
the DW velocity and shape, as well as for the depinning 
and critical (Walker breakdown) current densities. 

Ferromagnetic nanowires are natural candidates for 
building ultrafast memory and logic devices that rely 
on nanoscale current-induced DW motion. While quite 
high DW velocities (« 100 m/s) are possible in permalloy 
(NiFe) nanowires these setups often suffer from lim- 
ited reproducibility, strong DW pinning, and low critical 
currents. Following earlier proposals 10, 11|, recent ex- 
periments [13] realized Co nanowires in an AlO^; /Co/Ft 
trilayer structure, where the structural inversion asym- 
metry causes an interfacial electric field and thus a 
Rashba spin-orbit interaction (SOI) This electronic 
SOI is strong and allows to largely circumvent the above 
problems [13, [13, [l3| . The observed DW velocities of up 
to 400 m/s (and other interesting features, e.g., DW mo- 
tion against the current direction) were attributed [l2j 
to a conspiracy of theS-term and a field-like adiabatic 
Rashba spin torque [ifl [H [SS], 

denoted Ti below. 



back to the electronic bandstructure, it depends on the 
current and can switch the magnetization. These ex- 
citing experimental observations and their techn olog ical 
promise have triggered further theoretical work 2F25j 
that has drawn a more complex picture, involving also a 
Sloncezwski-type nonadiabatic Rashba spin torque (de- 
noted T2 below) due to the interplay of SOI and spin 
relaxation. 

Given the complexity of this problem, we here aim 
at understanding current-induced DW motion in a 
ferromagnetic Rashba nanowire for the simpler one- 
dimensional (ID) limit. This limit allows for the analyt- 
ical calculation of the full current-induced spin torque. 
Previous experiments [l2j have used Co nanowires of di- 
ameter K, 500 nm, much thicker than the few-channel 
ID case studied here, but future experiments could ap- 
proach this ultrathin-wire limit. By numerical solution 
of the Landau-Lifshitz-Gilbert (LLC) equation, with the 
nanowire initially containing a Bloch DW, we predict ul- 
trafast DW velocities in current-pulsed setups, plus a rich 
variety of interesting other features. 

Model. — We study an ultrathin (ID) ferromagnetic 
Rashba nanowire along the x-direction, where the dy- 
namics of the magnetization, — Msn(x,i), with unit vec- 
tor n and saturation magnetization Mg, is governed by 
the LLC equation. 



dtn = — n X Hcff + an x dtn + T. 
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The effective field Hoff generates magnetic texture in the 
absence of a current, and the Gilbert damping parame- 
ter a depends on intrinsic material properties. The spin 
torque T{x,t) encapsulates all current-induced contribu- 
tions due to the exchange interaction of n with itinerant 
spin-polarized electrons. The very same electrons, how- 
ever, also experience the SOI. To calculate T, we start 
from the ID Sugawara representation [2^ for the Hamil- 
tonian i?ei describing the itinerant electrons, which is 
formulated in terms of the chiral spin current operators 



While Ti does not involve spin transfer and can be traced 
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with ID fcrmion operators Cp^a describing left (p = L = 
+) and right (p = i? — — ) movers for spin projection a. 
Here, er is the vector of spin Pauh matrices, : • • • : denotes 
normal ordering, we use units where ft = 1, and the spin 
density is s(a;) = 3l + Jr- Similarly, the chiral particle 
currents, = '4>,<^'^p,<^^ with expectation value jp, 
respectively, determine the density pc = Jl + jn and the 
charge current density Ic oc (j'l — jn) flowing through 
the nanowire. It is well established [l^l that without 
SOI, the universal low-energy spin Hamiltonian is Hq = 
{v/2) J dx : Jp- Jp :, where the dispersion relation has 
been linearized around the Fermi momentum (kp) and v 
is the spin velocity. The Rashba SOI corresponds to the 
single- particle contribution aj^{aykx—<Txky), where k^ ~ 
—pkp and the channel- mixing term oc ky is negligible 
for ultrathin nanowires (27| . Adding the sd exchange 
interaction (exchange coupling Agd), we arrive at 



Ho 



p=± 



dx Jp • (n - paiiey) (3) 



with the dimensionless SOI coupling an = 2ajikp / lS.sd- 
The current-induced spin torque entering Eq. ([l} is then 
given by 



T{x,t) 



'Asdn X (Jl + Jr) 



(4) 



where the average is over the electronic degrees of free- 
dom. To obtain T, we set up the equations of motion for 
the spin currents. Using their Kac-Moody algebra (26j . 
we find {p ~ L, R) 

dtJp - pvdx3p = - Asd Jp X (n - paRSy) (5) 

- y: /3('')A.d(jp-jr), 

u—i^ii 

where spin relaxation has been included phenomenologi- 
cally within the relaxation-time approximation. We iden- 
tify two competing relaxation mechanisms characterized 
by different stationary spin current configurations, Jp ^ 
and Jp , and different rates governing the relaxation into 
them, P^^^Asd and /3(")Asd, respectively; for simplicity, 
we assume isotropic (p-independent) rates. Relaxation 
processes of type (i) try to establish a stationary distri- 
bution through boundary conditions, where the nanowire 
is contacted by ferromagnetic leads that inject a spin- 
polarized current with polarization factor < P < 1 
[2^ . The stationary spin current is then given by 



(6) 



Processes of type (ii) instead describe spin relaxation to 
the stationary solution of Eq. ([S]), 



Jf (. 



(7) 



Nonadiabatic spin torque with spin-orbit coupling. — 
Assuming slow space-time variations of n{x,t), we now 
proceed by solving Eq. ([5l through an iterative (space- 
time) gradient expansion [5| . The zeroth-order result (no 
terms involving space-time derivatives) already contains 
the leading terms of the two Rashba spin torques men- 
tioned above, T'^°\x,t) = Ti + T2. Indeed, with the 
spin-polarized current density Is ~ PIc/{2e), some alge- 
bra along the lines of Ref. [5| yields 
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We here quote only the result to lowest nonvanishing 
order in or, but complete expressions are given in the 
supplementary material [s^l- We stress that both spin 
torques are already present without magnetic texture 
once Or Their structure agrees with previous re- 

sults [2ll425| . For the first-order spin torque, we find 



-(1) 



J2 {Aud^n + B^xi X d^n) - H)^>n x ey,(10) 



-fffl = ^ [KdyUy + h'^{n X dyXl) ■ Gy] , 

y—x^t 

with explicit expressions for A^.t, Bx.t, h^^t and h'^ j given 
in the supplementary material [30|. Additional higher- 
order contributions in the gradient expansion Q will not 
be taken into account here. It is useful to distinguish two 
families of contributions to T^^^ First, we have terms 
oc dx^t'o. and oc n x dx,ti^, which already appear without 
SOI. Indeed, for aR ^ 0, Eq. (jTU]) contains the coefficients 



A, 



(11) 



We recognize the adiabatic spin torque (the term) as 
well as the "/3-tcrm" (the term) [sj. Note that the 
/3-term is caused by relaxation mechanism (i). Includ- 
ing now also the Rashba coupling aR^ we encounter a 
renormalization of A^ and B^ which, however, does not 
cause qualitative changes in the DW motion. Similarly, 
the time-dependent terms [At and Bt) mainly renormal- 
ize the Gilbert damping parameter a. The second family 
of spin torque contributions consists of terms oc n x e^, 
corresponding to an effective local magnetic field. 
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with Ti = Tin X Gy specified in Eq. 



We refer to the 
field in Eq. (fT2)l as "nonequilibrium Rashba field" below. 
When expanding the full result [s^l for Hr to second 
order in the SOI. wc obtain 



Hr 
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Besides the well-known Rashba field oc urIs due to 

Ti [13, [Til, [iB-i^, we find the additional contribution 
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FIG. 1. (Color online) Velocity of the DW center vs propa- 
gation time for (a) different or and (b) different 7c [A/m^]. 
An increase of either or or Ic increases the Rashba field Hr. 
At the critical value H^^, we find a current-induced Walker 
breakdown and the velocity starts to oscillate. Further in- 
crease of Or or Ic increases the oscillation frequency. 
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FIG. 2. (Color online) Velocity of the DW center vs propaga- 
tion time for five current pulses of different length tp. After 
the pulse, the DW drifts for a certain time in either positive 
or negative direction (depending on the oscillation phase at 
the pulse end). The dashed curve shows the corresponding 
result for a steady-state current. 
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oc aj^lg- This term dominates for (3^ 
thus leads to a strong noncquilibrium Rashba field even 
for modest Rashba coupling a^. While nonadiabatic 
torque contributions [T2 in Eq. © and the term in 
Eq. ([TT])] appear due to relaxation mechanism (i), H^^ 
is caused by mechanism (ii). 

Numerical simulation. — To obtain explicit results for 
the DW motion, we have performed numerical simula- 
tions [SlJ. We study a Bloch-y DW created by the effec- 
tive magnetic field Hoff = Jd^n + KuyBy — K±nxex in 
Eq. ([T]), where J is an exchange coupling between local- 
ized moments; K and K± are anisotropy constants. Mea- 
suring length (time) in units of xq — 0.5 nm {to = 1 fs), 
we adopt the following parameter values: J = 5.2, 
K = 0.185, K± = 0.008, a = = = 0.06. 

This choice corresponds to Ta/CoFeB/MgO, a material 
with comparably strong SOI as in Co nanowires and well 
characterized parameter values [l7j . Furthermore, unless 
stated otherwise, we set Agd/w = 0.2, = 0.2, P = 1, 
and Is = 0.06. The latter value implies the charge cur- 
rent density /c = 1 A/m^. In addition, the density pc 
and the spin velocity v appear in the full expressions for 
T [s^l, where we take v = 2x 10'''m/s and pc ~ Id v. We 
have also repeated the calculation for other parameter 
sets to verify that the results shown below are generic. 

Steady- state current. — Numerical simulation of Eq. ([T]) 
with a time-independent current density Ic confirms that 
the DW motion is strongly influenced by the Rashba SOI, 
both concerning the instantaneous (momentary) motion 
and the asymptotic long-time behavior. Let us first dis- 
cuss the time evolution of the momentary velocity, see 
Fig. [TJ With increasing SOI strength a^, we observe a 
transition from purely translatory motion, with (asymp- 
totically) constant DW velocity, to a regime with su- 
perimposed oscillations in the time-dependent velocity. 
Our numerical analysis reveals that the oscillations stem 
mainly from the Rashba field Hr- In fact, the DW mag- 



netization is found to preccss around the j/-axis, which is 
the direction of the Rashba field. This behavior strongly 
resembles a field-'mdnccd Walker breakdown, where the 
oscillations appear once Hn exceeds a certain critical field 
^WB^ As is known from the field- induced case, this 
critical field depends on the perpendicular anisotropy, 
i/WB K^^. Within the regime Hr > iJ^^^, we 
observe that the oscillation amplitude increases with K± 
but remains independent of all other system parameters. 
Moreover, the oscillation frequency is proportional to the 
value Hr of the noncquilibrium Rashba field. Since Hr is 
affected by the current (not only by ur), the DW motion 
can be effectively controlled by a spin-polarized current 
Is, see Fig. [TJb). The impact of all other parameters 
on the DW motion is less pronounced and not discussed 
here. 

DW response to current pulses. — Experiments are of- 
ten carried out with current pulses instead of steady-state 
currents l3, 1^ ■ For rectangular current pulses of width 
tp, we find that current-induced DW motion closely re- 
sembles the behavior in the field-driven case above the 
Walker breakdown Here, the DW does not immedi- 
ately stop at the end of the pulse, but instead drifts for 
a certain time with nearly constant velocity, see Fig. [2] 
Interestingly, the drift direction depends on the actual 
phase of the oscillation at the end of the pulse. There- 
fore, even if the DW initially moves forward due to the 
current pulse, it may still end up in a backward position 
relative to its starting point. For short current pulses, 
the drifting mechanism can completely dominate the to- 
tal DW displacement. 

Global average DW velocity. — In the current-pulsed 
case, a straightforward way to determine the average DW 
velocity is to measure the initial and final DW position 
before and directly after the pulse, and then to divide 



this distance by the pulse width tp 



In experiments, 



however, it is often difficult to read out the DW posi- 
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FIG. 3. (Color online) (a) Momentary velocity of the DW 
center vs propagation time for a steady-state current, (b) 
Average DW velocity {vdw) = XDw{t —J- oo) /tp vs current 
pulse width tp with Ic [A/ra^]. Depending on the momentary 
velocity at the pulse end, the DW drifts in different directions. 
Especially for short pulses, this may strongly affect {vdw}- 




FIG. 4. (Color online) Global average velocity of the DW 
center (vdw) vs current density Ic for fixed pulse width tp. 
Since Ic determines the oscillation frequency, the pulse can 
end up in different oscillation states, causing positive or neg- 
ative drift. This can strongly change (vdw) and causes DW 
motion against the current direction for small Ic- 



tion right at the end of the pulse, and usually one ob- 
tains the final position only somewhat later. The ex- 
tracted average velocity thus coincides with the real one 
only if the domain wall stops instantly when the pulse 
ends. Figure [2] demonstrates that such assumptions are 
not valid in general: The final position, XDw{t ~^ oo), 
typically deviates strongly from the location right after 
the pulse, XDw{tp)- In the following, we study the exper- 
imentally more accessible "global" average DW velocity, 
(vdw) = XDw{t ^ oo)/tp. A major issue determin- 
ing {vdw) is the oscillation phase reached by the DW 
at time t = tp. With changing tp, the oscillatory DW 
magnetization ends up in different configurations asso- 
ciated with different drift directions and final positions, 
XDwit oo) = XDwitp) + a^drift- WMlc XDw{tp) ^ tp 
is determined by the pulse width, Xdrift can have ei- 
ther sign. For very long pulses {tp — >• oo), we have 
|a^DVK(^p)| ^ |a;drift|, and the conventional average ve- 
locity, x{tp)/tp, coincides with {vdw)- For short pulses, 
however, they may differ considerably. We find that the 
tp-dependence of the global average velocity {vdw) ex- 
hibits a sawtooth-like behavior, see Fig. [3l 

Tuning the pulse width tp is not the only way to con- 
trol the motion of the DW. For fixed tp, we can also 
change the oscillation frequency of the DW motion by 
changing the current density Ic, see Fig. ID We again 
find a sawtooth-like dependence of {vdw) on Ic, which 
roughly oscillates around the velocity found for an = 0. 
For small values of Ic, the DW even moves against the 
current fiow. Due to the strong /c-dependencc of (w_dh/), 
small changes in the current can here lead to completely 
different types of DW motion. 

Conclusions- — In this paper, we have shown that 



Rashba spin-orbit interactions can qualitatively affect 
the current-driven domain wall (DW) motion in ferro- 
magnetic nanowires. Their main influence is encoded 
in a nonequilibrium Rashba field Hr, which is respon- 
sible for magnetization precession and thereby for qual- 
itative changes in the domain wall dynamics. Remark- 
ably, we find that the conspiracy of spin relaxation and 
spin-orbit coupling is able to generate a dominant con- 
tribution to Hn even for weak Rashba couplings. As 
a consequence, above a critical current-induced Rashba 
field Hr > H^^, the DW dynamics is predicted to ex- 
hibit oscillatory features (similar to the field-driven case 
above the Walker breakdown) and a nontrivial depen- 
dence of the DW velocity on the current density. In fact, 
the DW motion can turn from a purely translationary 
into an oscillatory dynamics. When using pulsed currents 
with variable pulse width, even the direction of the DW 
motion can be controlled. All predicted effects should be 
observable with present-day experimental tcchnic^ues. 
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SUPPLEMENTARY MATERIAL 

The total spin torque can be written in terms of orders 
of derivatives of the magnetization n as 

T=^tW, (14) 

k 

where T^'^^ contains only terms proportional to the kth 
derivative of n in time and/or space. We find for the 
zeroth order spin torque 

T(0)=Ti+T2, (15) 



with 
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-aRny{a% - 1) ^ n x e^^ , 



T2 =/3^''^^ [-^.(1 + 4) + '^Psvanny] n x (n x e,,) , 
N ={l + alf - Aalnl , 

with ps = Ppc/{'2,e) and Ig = PIc/{2e). Ig = PIc/{2e), polarization P, charge density pc, charge current density Ic, 
elementary charge e and the inverse relaxation times (3^'^^ = (AsdT^"')""'^. The index ly ~ denotes relaxation due 
to the leads (i) or relaxation into the stationary state of the isolated wire (ii). 



Additionally we get for the first order contribution rather lengthy expressions as 



with 



and 
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